Two-Frequency Jahn- Teller Systems in Circuit QED 
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We investigate the simulation of Jahn- Teller models with two non-degenerate vibrational modes 
using a circuit QED architecture. Typical Jahn- Teller systems are anisotropic and require at least a 
two-frequency description. The proposed simulator consists of two superconducting lumped-element 
resonators interacting with a common flux qubit in the ultrastrong coupling regime. We translate 
the circuit QED model of the system to a two-frequency Jahn- Teller Hamiltonian and calculate its 
energy eigenvalues and the emission spectrum of the cavities. It is shown that the system can be 
systematically tuned to an effective single mode Hamiltonian from the two-mode model by varying 
the coupling strength between the resonators. The flexibility in manipulating the parameters of the 
circuit QED simulator permits isolating the effective single frequency and pure two- frequency effects 
in the spectral response of Jahn- Teller systems. 

PACS numbers: 42.50.Pq, 71.70.Ej,85.25.-j 
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I. INTRODUCTION 

Simulating complex physical phenomena using systems 
that offer precise control of physical interactions, such as 
ultracold atoms [l[ , Bose- Einstein condensates @, H| and 
trapped ions [3| has attracted much attention in the last 
decade. It has been recently shown that cavity QED 
systems can be utilized for the same purpose, in partic- 
ular to simulate certain gauge potentials, the anomalous 
Hall effect, and the Dirac equation [f|. The potential use 
of cavity QED systems to simulate such physical mod- 
els relies on the successful simulation of Jahn- Teller (JT) 
interactions which require atom-photon ultrastrong 
coupling conditions [9[ . 

JT models describe the interaction of localized elec- 
tronic states with vibrational (phonon) modes in crystals 
or in molecules 10]. Cavity QED systems have been al- 
ready proposed to simulate single mode JT models [TTj . 
On the other hand, many practical systems need a de- 
scription in terms of multi-mode JT interactions [TH-EtJ . 
We address the question of how to generalize the restric- 
tive single mode simulation of JT systems to two-mode 
JT interactions within the circuit QED context. 

Circuit QED 18| offers the possibility to operate in 
the ultrastrong coupling regime [l9l - |24j for efficient JT 
and related spin-boson or Dicke model simulations. Pho- 
tonic waveguide arrays are alternatively proposed [25| for 
reaching the deep ultrastrong coupling regime (DSC) [26[ 
of JT type light-matter interaction. 

Our idea is to consider a system consisting of a two- 
level atom simultaneously interacting with two cavities 
coupled to each other, rather than with a degenerate two- 
mode cavity, which was considered for th e cy lindrically 
symmetric E x e JT model in cavity QED [ll[. In terms 
of the normal modes of the coupled cavities, our system 
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allows simulating a two-frequency (two non-degenerate 
vibrational nor mal modes) E X (ft + ft) JT model [ifjj . 
The normal modes of the two coupled cavities consist of 
a high frequency component and a low frequency compo- 
nent. The coupling strength between the cavities can be 
utilized to alter the frequency ratio of the modes to sim- 
ulate different frequency ratios encountered in different 
JT impurities in solids [27] . In addition to more realistic 
simulations of JT systems, establishing a link between 
multi-mode JT models and coupled circuit QED systems 
could enable exploring many body physics such as quan- 
tum chaos [13, EU , quantum phase transitions [l6[ , and 
quantum entanglement in JT systems [24| [28| by using 
coupled cavity arrays. 

We consider a coupled circuit QED [TJ, [29rl33l ] system 
in the ultrastrong coupling regime as a feasible platform 
on which to realize our idea. The system consists of two 
coupled lumped-element LC resonators interacting with 
a two- level artificial atom, a superconducting flux qubit. 
In the ultrastrong coupling regime the rotating- wave ap- 
proximation is not valid, so that the qubit-resonator cou- 
plin g is of JT type rather than Jaynes-Cummings type 
[221 |24| . allowing strongly coupled multi- frequency JT 
systems to be simulated. The switchable ultrastrong cou- 
pling architecture can also be applied if a tunable cou- 
pling strength between the resonators and the flux qubit 
is desired [34]]. 

Typical treatments of strongly coupled multi-mode JT 
systems in chemistry or in condensed matter physics uti- 
lize a cluster model 35], or use an effective single mode 
model where most of the JT interaction energy is concen- 
trated predominantly upon a single effective vibrational 
frequency with a negligible spread (narrow range of fre- 
quencies) [36[ . These methods are especially used for in- 
terpreting effects associated with low energy states, such 
as those seen in low temperature optical absorption [37| . 

When the frequency difference between the two modes 
is large, the situation is analogous to the case of optical 
and acoustic phonons in solids for which perturbative cor- 
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rections become significant. We show that the frequency 
separation of the modes over which the JT interaction is 
distributed can be tuned with the coupling strength be- 
tween the resonators. Our coupled circuit QED proposal 
allows for systematic simulation of effective single mode 
and pure two-mode effects as well as transitions between 
these regimes. 

This manuscript is organized as follows. In Sec. [TTI we 
introduce the two-mode JT model and its implementa- 
tion in a circuit QED context. The effective single mode 
treatment is described in Sec. IIIII In Sec. IIV1 the ex- 
perimental implementation is laid out. The results and 
discussions are presented in Sec. [V] Finally we give con- 
clusions in Sec. I VII 

II. TWO-FREQUENCY JT SYSTEMS IN 
CAVITY AND CIRCUIT QED 

When two electronic levels are coupled to vibrations 
of ions or atoms in solids and molecules the general form 
of the interaction can be written as Hjt = f(Q) ■ <?, 
where f(Q) is a vector valued function of vibrational co- 
ordinates while a is the vector of the Pauli spin matri- 
ces. Such interactions are in general called Jahn- Teller 
interactions [lfj. In this paper, we focus on a particular 
one of the form (Q± x + Q2x)&x, which is known as the 
E x {fix + fa) JT (or Herzberg- Teller) model. Our aim 
is to generalize the recently discussed simulation of sin- 
gle mode E x e JT model in cavity QED [III 1 which 
is of the form Q x a x + Q y a y , to the two-frequency case. 
Our choice of E X {(3x + $2) is the simplest possible two- 
frequency JT model. The implementation of Ex (/3i+/?2) 
model allows for simulating realistic crystals that exhibit 
spatial anisotropy. The single-boson E x (3 model is for- 
mally equivalent to the Dicke model and signatures of 
quantum chaos have been discussed in the E x {j3\ + P2) 
model recently (l6| . 

The Hamiltonian corresponding to the multi-frequency 
JT interaction between a single impurity ion and many 
vibrational degrees of freedom of the host lattice (or 
molecule) is expressed as 

H = H p h + Hjt , (1) 

where H p h describes the free Hamiltonian of the phonon 
modes at frequencies u>i 

H ph = y^i(alai), (2) 

i 

where ai(a|) are the annihilation (creation) operators of 
the phonons. The multi-mode JT interaction describes 
the coupling of the single ion to the vibrational modes 

Hjt = £)fiw i fc ( (a i + aJ)V. ( 3 ) 

i 

Here ki are the dimensionless scaling factors of the JT 
coupling coefficients and V is an operator that depends 
on the electronic degrees of freedom of the impurity ion. 



We wish to simulate this multi-mode JT interaction 
using a coupled two-resonator circuit QED system. Nor- 
mal modes of the coupled microwave photons play the 
role of phonons, while a flux qubit plays the role of the 
impurity ion. The interaction of the flux qubit in the two- 
resonator circuit QED system mimics the local (short- 
range) interaction of the ion-phonon coupling. On the 
other hand, there is an additional non-local (long range) 
coupling between the resonator modes, describing hop- 
ping of photons between the resonators in the circuit 
QED system which mimics the coupling between the vi- 
brational phonons. The coupled resonator model can be 
written as {h = 1) 

H = H q + H c + H qc + H cc , (4) 

where 

Q 

H q = ^<r g , ( 5 ) 

H c = riidjai + f^a^d^, (6) 
H qc = + ax) + A 2 (a| + a 2 )}o x , (7) 

H cc = J{dx*d2 + d 2 'dx), (8) 

where a z , and a x — 0+ + er_ are the Pauli spin operators 
describing the qubit degrees of freedom with Q being the 
qubit transition frequency, Qx.2 the resonance frequen- 
cies of the cavity modes, J the hopping rate of microwave 
photons between the resonators and &x,i {&\ 2) the anni- 
hilation (creation) operators for the cavity photons. We 
want to simulate the two-frequency JT model in Eq. [1] 
with the two-resonator circuit QED model in Eq. @] For 
that aim it is necessary to be able to transform one model 
to the other and show that they are identical for a certain 
set of model parameters. In the next section we examine 
the transformation between these two Hamiltonians. We 
apply the so-called effective single privileged mode trans- 
formation that has been developed for multi-frequency 
JT systems \^ that allows for systematic analysis of 
pure single and multi-frequency effects. 

III. EFFECTIVE SINGLE-MODE JT SYSTEM 
IN TWO-RESONATOR CIRCUIT QED 

We now employ the effective single mode treatment 
[36| for the two-frequency JT model obtained within the 
two-resonator circuit QED context. For that aim we look 
for a particular linear superposition the of normal modes 
of the coupled resonators, 

dj = ^ Ajk&k, (9) 

k 

where Aik are the elements of a real orthogonal matrix 
to be determined, such that most of the JT energy is 
concentrated over a privileged mode among the set of 
new bosonic modes dj. 
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Without loss of generality, we choose the privileged 
mode as di for which the single-mode model JT system 
can be written in terms of an effective frequency u e f f and 
an effective JT coupling k e f / , scaled by u) e f / , so that 

H eff = 7^ + Ueff[a{&i + fc e //(di + a\)V]. (10) 

The relevant elements of the transformation, An and 
A\2, are determined by maximizing k^ffOJeff- This is the 
amount by which the minimum of the potential energy 
of the system is lowered, under adiabatic approximation, 
due to the interaction of the rest of the system with such a 
single-mode [3l| . This gives An : A\2 = k\ : fca, subject 
to normalization conditions. Direct substitution of Eq. [5] 
into Eq. Q] yields 



k 2 



K eff 

k 2 . 



(11) 
(12) 



For a 2 x 2 orthogonal matrix A determining the first 
row of elements fixes the remaining two by orthonormal- 
ity conditions such that 
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fci k 2 
k e ff V -ki 



(13) 



which is taken in consistency with Ref . [36| . Such trans- 
formations are common in Morris-Shore bright and dark 
state transformations [38| . For the system analyzed here, 
on the contrary, there is no perfect decoupling of cither 
mode from the dynamics of the rest of the system, though 
for certain parameter regimes the modes di and ai be- 
come approximately decoupled. The total transformed 
Hamiltonian can be written as 



H S y S h, 



eff 



H'ph 



Hi; 



where 



with 



H ' P h 



LOlkn 



UJ 2 kj 



(14) 
(15) 

(16) 
(17) 



is the free Hamiltonian of the disadvantaged effective 
mode. Interaction of this mode with the rest of the sys- 
tem is described by 

Hint = c 2 [(a\a 2 + ax&l) + k eff (a 2 + a\)V]. (18) 

Here the strength of the coupling between the privileged 
and the disadvantaged modes is characterized by the pa- 
rameter 



C2 



k x k 2 

K eff 



(19) 



where A = u\ — co 2 is the frequency difference between the 
vibration modes in the two-frequency JT model. If the 
vibrational modes are degenerate the model is exactly 
equivalent to the case of an effective single frequency. 
The coupling between the JT vibration modes can also 



be expressed as c 2 



j 2 — cj , where 
Lo r {k\+u:^kl 



(20) 



This allows for interpreting c 2 as the mean square width 
of that distribution with a mean ZU ~ uj e f f [361 ] . 

The effect of Hi nt on the effective single-mode model 
can be examined perturbatively provided that the fre- 
quency spread (A or c 2 ) is not too large. Perturbative 
effects will only be significant on the JT ground state 
only starting from the second order. To see this it is 
convenient to introduce a new set of operators [36[ 

fj = a 1 +k eff V, (21) 

for which we can re-express the effective model as 



H P ff = i 



l eff ^^effV'V-^effkl ff V 2 . (22) 

The last term is proportional to the unit matrix, i) obeys 
bosonic commutation rules even though it contains Pauli 
spin operators. This operator can be identified as a 
"bright" qubit-polariton quasi-particle. The first term 
in the effective Hamiltonain is then a harmonic oscilla- 
tor contribution due to the free energy of these quasi- 
particles. The interaction Hamiltonian becomes 



Hint = c 2 {r)6t\ 



(23) 



representing the coupling of a bright qubit-polariton to a 
"dark" effective mode a 2 . Using commutation relations 
between H e ff and a.\, the necessary matrix elements of 
f) for the perturbation analysis in the representation of 
eigenstates of H e ff can be determined [36j. The simple 
relation [-ff e //, di] = —w e fff] shows that 77 has no di- 
agonal elements. This means that the perturbative cor- 
rections can only be of significance starting from second 
order, at least for the JT ground state and for low energy 
lying states. If the frequency spread c 2 is not negligible 
compared to the Jahn- Teller coupling fc e // the descrip- 
tion of the system using a privileged single mode is not 
possible since the JT energy is spread among the two 
modes di, a 2 . 

We are now in a position to relate the two-resonator 
circuit QED model of Eq. 0] to the two-frequency JT 
model of Eq. [TJ We find that the parameters of the 
Hamiltonians are related as follows 

W e// = ^1, U =0,2 W e ffk e ff = \l, 

c 2 k eff = A 2 , c 2 = J. (24) 
The relations require that a condition of the form 



n 1 = ±j 



(25) 



should be satisfied among the parameters of the circuit 
QED Hamiltonian. 
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IV. EXPERIMENTAL IMPLEMENTATION 

The system we consider to implement the Hamiltonian 
in Eq. @] consists of two lumped-element LC resonators 
capacitively coupled to each other and a flux qubit simul- 
taneously coupled to each resonator. A schematic of the 
circuit can be seen in Fig. [TJ The resonator-resonator 



a) 




b) 



FIG. 1: a) (Color online) Schematic of a flux qubit gal- 
vanically coupled to two LC resonators through coupling 
inductances L C1 2 . Each resonator has resonant frequency 

£^1,2 = ((£1,2 +L C1 2 )Ci,2) -1 ^ 2 . The resonators are coupled to 
each other through the coupling capacitor C c . b) Implemen- 
tation of the circuit using interdigitate finger capacitors. Feed 
lines individually coupled to each resonator can be employed 
to probe the internal photon state of the resonators. 



interaction J is determined by the coupling capacitor 
C c between the two LC resonators when they are in 

the ground state J ~ C c V IinSl V Tm 



a 



where 



4CiC 2 ' 

Ci } 2 is the total capacitance of the uncoupled resonators. 
Using typical sample parameters pol . [23 |. the coupling 
strength between the two resonators can be made very 
large, up to a considerable fraction of the frequency 
of each resonator. The spurious inductive resonator- 
resonator coupling could be minimized with an appro- 
priate resonator geometry, if necessary. The flux qubi t 
can also be made of large enough size (as the ones in |4fJ| ) 
so as to increase the distance between the resonators and 
reduce the mutual resonator-resonator inductance. 

The coupling energy A in a flux qubit-resonator system 
can become a large fraction of the energy of the resonator 
if the qubit is galvanically attached to the resonator, as 
already demonstrated experimentally in Refs. [ljj [20| . 
For a qubit either sharing a long section of its inductance 
[20l | or a Josephson junction [l9|, [22j with a resonator, 
coupling energies A rs cdj.,2 are within reach experimen- 
tally. 

From the analysis of the sec. IIII1 in order to study 
the privileged mode regime the coupling term in Eq. [TU] 
needs to be larger than the coupling term in Eq. 1181 
This implies, according to the relations in Eq. [5H that 



Ai > A2,J. Therefore in the experiment the flux 
qubit has to be ultrastrongly coupled to one resonator 
and strongly coupled to the other resonator, while the 
resonator-resonator coupling must be close to the qubit- 
resonator strong coupling. These designed coupling en- 
ergies will determine the privileged mode. By detec ting 
the photon state of each resonator using feed lines [4J| 
permits exploring the spectral properties of the complete 
system. 



V. RESULTS 

For the sake of simplicity, we choose fci = k 2 = k so 
that k e ff — \/2k. The relations in Eq. |2"41 reduce to 
fix = fi 2 = ^ c = (wi + w 2 )/2, Ai = (wi + u 2 )k/V2, A 2 = 
Ak/y/2, and J = c 2 = A/2. Our choice requires the 
resonators in the circuit QED system to be degenerate. 
We further assume resonance condition £1 — fl c . The 
circuit QED Hamiltonian then becomes 



H 



a[ai + a 2 a 2 + -a z + — (a[a 2 



■ a\dt\) 



A 



+ k e ff[(a[ + — {a\ + a 2 )](Tx- 



(26) 



We use dimensionless energy and time, respectively 
scaled by Ml c and l/Q c > but do not change our notation 
for scaled variables. Our model is then a two-parameter 
(fc,A) theory where A is in units of fl c . 

The resonators are degenerate but the system still sim- 
ulates the two-frequency JT model. The coupling coeffi- 
cient between the resonators J determines the frequency 
ratio of the two vibration modes in the corresponding 
two- frequency JT Hamiltonian of Eq. [T]to be simulated, 
which becomes 



H 



1 



wia[ai + w 2 a 2 a 2 + -a z + 
+ k[uii(a\ + fix) + uj 2 {al + a 2 )](J x , (27) 



where the frequency ratio is determined by 



uji _ l + A/2 
~ 2 ~ 1 - A/2' 



(28) 



Some cases of interest are the 2:1 frequency ratio of the 
two phonon modes in CgHg 1 * 1 , and the frequency ratio 



27J. As 



3 : 1 of the two phonon modes of Fe 2+ in ZnS 
explained in sec. Mil the transformation or equivalence of 
these two Hamiltonians is exact. For A (or the frequency 
spread c 2 ) small compared to the uncoupled eigenfre- 
quencies of the qubit and the resonators, we obtain a 
faithful representation of the two-frequency JT model in 
terms of the privileged mode. Corrections appear only 
as second order perturbations. 

We first examine the eigenenergies of the Hamiltonian 
in Eq.[5U The lowest five eigenvalues are shown in Fig. [21 
Here the Fock space dimensions for each resonator mode 
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FIG. 2: (Color online) Dependence of the lowest five eigenen- 
ergies of the two-frequency Jahn- Teller model on the mode 
coupling parameter A. (a) For the qubit-cavity coupling scal- 
ing parameter k — 0.1/\/2, the splitting of the first energy 
level into three at A = is a pure Rabi splitting. The privi- 
leged effective single mode regime is valid over a narrow band 
A < 0.1. (b) At k = 1/V2 the energy-level transitions 
become flat over a broader band |A| < 1, deeper into the 
ultrastrong coupling regime, (c) Energy bands of the two- 
frequency Jahn- Teller model depending on the mode coupling 
parameter A and the qubit-cavity coupling scaling parameter 
k. All the quantities plotted are dimensionless as explained 
in the text. 



is fixed to 2 so that we consider up to two-photon man- 
ifolds. We examined the influence of dimensions of the 



Fock space on our results in the case of spectrum cal- 
culations and found them to be sufficiently robust. The 
strong coupling case with k = 0.1/ \/2 (or just the be- 
ginning of the ultrastrong coupling regime) is considered 
in Fig. 2(a) When A — there is only pure Rabi split- 
ting as can be seen in the first excited level. When A 
increases, the coupling between the privileged and the 
disadvantaged modes increases. This polaritonic interac- 
tion of the modes causes further repelling of the Rabi- 
split levels. Single privileged effective mode description 
of the system would only be valid over a narrow band 
k ~ |A| < 0.1. The ultrastrong coupling regime with 
k = l/\/2 is considered in Fig. 2(b) An asymmetric 



Rabi splitting at A — can be seen in the first excited 
level. Here the effective single mode description is valid 
over a broader range |A| < 1. The dependence of the 
energy spectrum on the full range of k and A is shown 
in Fig. 2(c) The first band is tent-shaped, and for low k 



it varies sharply with A resulting in a narrow regime of 
the effective single mode description. As k reaches ultra- 
strong coupling conditions, the regime of effective single 
mode description becomes more robust against variations 
in A over a broader range. 

Solid state and molecular multi-frequency JT systems 
are usually investigated through their absorption spec- 
trum. The corresponding quantity in circuit QED is 
the transmission spectrum of the resonators. We con- 
sider the power spectrum of only one resonator, corre- 
sponding to the privileged mode. Deviations from sin- 
gle mode behavior in this spectrum would be identified 
as pure two-frequency effects. In order to calculate the 
power spectrum it is necessary to solve the quantum mas- 
ter equation for the ultrastrong coupling regime, which 
can only be rigorously formulated in the dressed state 
picture of coupled qubit-resonator, examined recently in 
refs. [H, [43[. Our purpose is to see quantitative changes 
in the spectrum at different frequency ratios of the two- 
frequency JT model. We assume that the usual Bloch- 
Redfield quantum master equation for circuit QED sys- 
tems in the Born-Markov approximation is applicable for 
and the equation is given as (h = 1) 



our purposes 



41 



^ = -i[H,p]+Cp, 
where the Liouvillian superoperator C is given by 
Cp = (1 + nth)KT>[aj]p + n th KD[a\]p 

+ ^V[a}p+^-V[a z }p, 



(29) 



(30) 



with nth being the average thermal photon number, 
which we take as n t h = 0.1 corresponding to 100 mK 
[4lT ] . The Lindblad type damping superoperators are de- 
noted by V. The cavity photon loss rate k is taken to be 
the same for both resonators. The qubit relaxation and 
dephasing rates are represented by 7 and 7^, respectively. 

The JT spectrum is determined by calculating the real 
part of the Fourier transform of the stationary two-time 
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first-order correlation function for the privileged mode 
&i, so that 



p(u) = / (d t 1 (t)d 1 (o)} ( 



(31) 



We use Python programming language with the QuTip 
package for the determination of the spectrum The 
decay parameters, scaled by f2 c , are taken to be k = 
0.001,7 = 0.001 and 70 = 0.01, while assuming Fock 
space of up to two photons for each resonator mode. 
When we consider higher Fock space dimensions (due to 
numerical constraints we examined up to 5 particle man- 
ifolds) we find that the spectrum is robust against the 
variations in the Fock space dimensions for the regimes 
of J < 0.5 we consider. At larger J values small changes 
in the spectral intensities are observed, but they are still 
negligible up to J ~ 1. For smaller decay rates and at 
even larger J values the spectrum becomes more sensitive 
to dimensions of the Fock space. 

Our results for different values of J are presented in 
Fig. [U Fig. 3(a) shows the spectrum when the two 
resonators are uncoupled, J = 0. For low k e ff, the 
spectrum shows typical asymmetric Rabi-split frequency 
peaks of the Jaynes-Cummings (JC) model around the 
degenerate frequency of the resonators. At larger values 
of k e f f ~ 1 the system is in the single-mode JT regime. 



Fig. |3(b) presents the effect of coupling the resonators 
with J = 1/2, which corresponds to the typical frequency 
ratio 3 : 1 of two phonon modes of Fe 2+ in ZnS. The po- 
laritonic splitting here shifts the Rabi-split peaks further 



away at low values of k, 



eff- 



For k e 



1 the system 



is in the two-mode JT regime. Fig. 3(c) shows similar 
features for J = 1/3, which corresponds to the typical 
value of frequency ratio 2 : 1 of the two phonon modes 
in CqKq^. In this case for k e ff > 1.5 the system enters 
the effective single mode JT regime. 

The general behavior of the spectrum with the 
resonator-resonator coupling J at a given qubit-resonator 
coupling k e ff is shown in Fig. 2J Fig. 4(a) shows the 
case when k is near the threshold of the ultrastrong cou- 
pling regime, k e ff = 0.1. The spectrum is mainly de- 
termined by the two-frequency character of the system, 
corresponding to the two-mode Jaynes-Cummings model. 
The privileged single mode description is limited to very 
small coupling strengths J = A/2 < 0.1. Beyond this 
point, normal mode splitting increases with J. 

When we consider the case of deeper ultrastrong cou- 
pling k e ff = 1, Fig. 4(b) reveals that it is easier to re- 
solve the single- mo de / two- mo de regimes as the single- 
mode regime is significantly enhanced, up to J = A/2 ~ 
0.4 < k e ff- This observation complies with our previ- 
ous arguments, based upon the energy levels of the sys- 
tem in Fig. [31 Beyond the single mode regime, Fig. |4(b)| 
shows that the higher frequency peak at uj ~ 1.4 disap- 
pears, while the lower frequency one at uj ~ 0.2 domi- 
nates. The spectrum exhibits additional peaks that grow 
in number and in amplitude in the two-frequency regime. 
These peaks are due to multi-photon processes that be- 






FIG. 3: (Color online) Resonator-qubit coupling k e ff depen- 
dence of the cavity emission spectrum for the privileged effec- 
tive mode in the two-frequency Jahn- Teller model for differ- 
ent values of the resonator-resonator coupling strength J. (a) 
J — corresponds to the pure single-mode regime. fe e // <C 1 
is the single-mode regime in the JC model. For k e ff > 1 
one enters the regime of single- mode JT model, (b) J = 1/2 
shows three regimes: k e ff <S 1 (2- mode JC model), k e ff > 1 
(2-mode JT model), and fc e // ^> J (effective single-mode JT 
model). This case simulates the frequency ratio of 3 : 1 of the 
two phonon modes in Fe 2+ in ZnS. (c) J = 1/3 simulates the 
frequency ratio of 2 : 1 of the two phonon modes in CeHe . 
In this case the effective single-mode regime is more clear for 
keff > 1.5. 
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(b) 

FIG. 4: (Color online) Variation of cavity emission spec- 
trum for the privileged effective mode in the two-frequency 
Jahn- Teller model with resonator-resonator coupling J at (a) 
k e ff — 0.1 and (b) fc e // = 1- 



cuit QED system. The proposed model consists of a flux 
qubit coupled to two resonators in the ultrastrong cou- 
pling regime. An exact transformation between the two- 
frequency JT Hamiltonian and the circuit QED Hamilto- 
nian has been established. The transformation permits 
describing the system in terms of an effective privileged 
single mode under certain conditions of the control pa- 
rameters of the circuit QED system. The effective dis- 
advantaged mode can be de-coupled from the privileged 
one in the ultrastrong coupling regime. The eigenenergy 
spectrum and power spectrum are calculated using ultra- 
strong circuit QED parameters, with specific attention to 
the present experimental restrictions. The tunability of 
the pure two-mode JT model and the effective privileged 
mode model is found to be feasible in the ultrastrong 
coupling circuit QED within the range of parameters in 
present experiments. 

Simulating and interpreting more complex JT sys- 
tems, such as vacancies in graphite or fullcrides C^~ , 
would require going beyond a two-mode description. 
Our analysis of the two-frequency JT model simulation 
can serve as a building block for further realizations of 
other classes of multi-mode JT systems, by consider- 
ing for example coupled multi-mode superconducting 
transmission line resonators and their interactions with 
flux qubits in the (ultra)strong coupling regime. Such 
extensions of the present work would allow examining 
rich geometric phase effects [l2| and designing synthetic 
gauge fields, as well as enhancing the comprehension of 
nonlinear JT dynamics of complex molecular systems. 



come more and more significant as one goes deeper into 
the strongly coupled JT model [45|. The higher energy 
resonance in the single mode regime turns out to be more 
susceptible to such multi-photon processes. The ampli- 
tude of this transition decreases and eventually vanishes 
in the two-frequency regime, while the lower energy reso- 
nance is more robust and does not decrease its amplitude 
significantly. These results suggest that one can monitor 
and analyze the transition between the effective single 
privileged mode and the pure two-frequency behavior of 
a JT system by tuning the circuit QED parameters into 
the DSC regime. 

VI. CONCLUSION 

In summary, we have presented a method to simulate 
a two-frequency JT model by using a two-resonator cir- 
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